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<V = (dv ~ 2) /2 of the muon where the role of input data needed in the evaluation 
is lowered in the interval between 1.2 and 3.0 GeV below charm threshold. This 



is achieved by decreasing the size of the weight function in the dispersion integral 
over the experimental ratio R(s) by subtracting a polynomial from the weight func- 
tion which mimics its energy dependence in that given energy interval. In order to 
compensate for this subtraction, the same polynomial weight integral is added again 
but is now evaluated on a circular contour in the complex plane using QCD and 
global duality. For the hadronic contribution to the shift in the anomalous mag- 
netic moment of the muon we obtain a^ ad (LO) = (701.3 ± 6.4) x 10~ 10 at leading 
order in the electromagnetic coupling. In addition, using the same procedure, we 
recalculate the next-to-leading contribution a)f d (NLO) = (-10.3 ±0.2) x 10~ 10 . In- 
cluding QED, electroweak, and light-by-light contribution, we obtain a value = 
(11 659 185.6 ± 6.4 had ± 3.5 LBL ± 0.4 QED+E w) x lO" 10 . 



1 Introduction 

The anomalous magnetic moment of the muon is one of a few physical parameters 
which can be determined with high precision and therefore can serve as precision test for 
the Standard Model of elementary particle physics. In this paper we join the attempts 
to determine the hadronic contribution of the anomalous magnetic moment of the muon 
with a better accuracy (for an overview over the status of calculations of different collabo- 
rations two years ago see e.g. Ref. PUI2])- While the dominant QED contribution a^ ED = 
(11658470.6±0.3) x KT 10 and the weak contribution a™ cak = (15.4±0.1±0.2) x KT 10 are 
known with very good accuracy (see e.g. Refs. [H1IUE] arid references therein), the main un- 
certainty is given by the hadronic contribution. As an example for the actual calculations of 
the leading order hadronic contribution we cite the value a^ ad = (683.1±5.9±2.0)xl0 -10 [§]. 

The calculation of the hadronic contribution to the anomalous magnetic moment of the 
muon mostly relies on experimental data in the e + e~ channel. In principle one can also 
make use of data from r-decay [7 . However, the inclusion of r-decay data introduces sys- 
tematic uncertainties originating from isospin symmetry breaking effects which are difficult 
to estimate fS|- We have therefore decided to only include e + e~ data in our analysis. 

There are different concepts for using the e + e~ data sets, either to its full extent [H] or 
in part by dividing the energy range into resonance regions close to the pair production 
thresholds where the experimental values are taken and those regions far from thresholds 
where perturbation theory is assumed to be valid. In this paper we present an approach 
which has been successfully applied for the determination of the running fine structure 
constant at the Z° boson resonance, ct(M§) [TU] . 

The subject of a precision determination of the anomalous magnetic moment became 
rather important again as a precision measurement for the positive muon could be accom- 
plished at the Brookhaven National Laboratory (BNL) jTT]. The value 
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cited in Ref. jTT] showed a deviation from the Standard Model prediction at that time, 

af 1 = (11 659 159.7 ±6.7) x 1(T 10 (2) 

as weighted average over the calculation results of different collaborations [2]. The devia- 
tion was 2.6 standard deviations which seemed to have opened the window for new physics. 
Many such suggestions were published, including concepts of supersymmetry, leptoquarks, 
lepton number violating models, technicolor models, string theory concepts, extra dimen- 
sions and so on. The discrepancy between measurement and the Standard Model prediction 
became smaller, though, when a sign error was discovered in the theoretical calculation of 
the light-by- light contribution Therefore, it is worth still to analyze the situation of 
the Standard Model prediction thoroughly, giving special emphasis on the error estimate, 
in order to compare it with the experimental world average 

a^P = (11 659 203 ±8) x 1(T 10 (3) 

which is dominated by the measurement of Ref. ^2] , as we will do in this paper. 



2 Theoretical background 

The leading hadronic contribution to the anomalous magnetic moment of the muon a^ ad 
can be extracted from the 0(a 2 ) vertex correction shown in Fig.^ The analytic expression 
for a^ ad from this diagram is a weighted dispersion integral over the imaginary part of the 
vector part of the hadronic vacuum polarisation, 

a had = 4o^ f 00 lmTl ^ {s) ^l ds (4) 
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with the QED kernel 
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Figure 1: 0(a 2 ) correction to the hadronic part of a M . The shaded bubble represents 
the current contributions to the vector part of the hadronic vacuum polarisation given in 

Refs. nang. 



As usual we make use of the kinematic variables 

(3(s) - 1 



x(s) 



/3(s) + Y V s 

ImIT had (s) is calculated from the hadronic current-current correlator, 
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Higher perturbative QCD corrections are well-known up to 0(a 2 ) with 0(m 12 / q 12 ) 
quark mass corrections which are supplemented by massless terms up to O(af) This 
perturbative part is denoted by U p (q 2 ). For brevity we only cite the first few terms of this 



expansion in a s and m 2 /q 2 
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with L = \n(/jr/q ). In our analysis we use the full 0(a s: m jq ) expression given in 
Ref. [T3|. The number of active flavours is denoted by n/ which changes according to the 
energy interval under consideration. For the zeroth order term in the m\ 2 Jq 2 expansion we 
have included the known higher order terms in a s , 
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with fco — 1, ki — 1.63982 and &2 = 6.37101. We have denoted the yet unknown con- 
stant term in the four-loop contribution by C3 which, however, does not contribute to our 
calculations since it has no absorption part. 

For the nonperturbative part the operator product expansion leads to 

T^Pf„2\ _ 1 ft i ^ a s\ / a s 

y2\2 



a \ 2 / a 
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36vr 2 (g 2 ) 2 V v ■ ■ / n j 2 43(g 

1 -(0.48 GeV) 8 - -^(0.3 GeV) 10 + O((g 2 )- 6 ) (10) 
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where Cp = 4/3, Ca = 3, Tp = 1/2 are 577(3) colour factors. We will use these results for 
the evaluation of the theoretical contributions. 



3 Experimental contributions 

The evaluation of the integral involves experimental data via the optical theorem which 
connects the imaginary part of the vacuum polarisation with the ratio 

R(s) = g( f f - ha ^°" S) = 12vrImn had ( S ). (11) 
cr(e + e — > ) 

Especially in the MeV and the low GeV energy range and in regions of resonances where 
perturbative QCD cannot be applied, data sets taken from experimental measurements are 
crucial. For our evaluation we have used the combined e + e~ data sets from Ref. [3 EH] 
and complement them in the dominant low energy range from 610 MeV to 961 MeV by 
recent two pion data By comparing with the new data for three and four pion 

decays f7[ E| one ensures that their influence for the sub GeV-range where the two pion 
data were recorded is neglible. 



All data sets are combined by weighting their relative contribution to the total experi- 
mental value with the corresponding statistic and systematic errors. 



4 Introduction of the method 

Global duality states that QCD can be used in weighted integrals over a spectral function 
if the spectral function is multiplied by polynomial functions. However, this does not work 
if the polynomial function is replaced by a singular function such as the weight K(s) /s in 
the present case. Nevertheless, local duality is expected to hold for large values of s far 
from resonances and threshold regions, i.e. Imll had (s) ~ Imil^ CD (s). 

In our approach we attempt to reduce the influence of experimental R(s) data in regions 
where the data has large uncertainties. Specifically, this holds for the interval between 1.2 
and 3.0 GeV. The essence of our method is to diminish the magnitude of the weight 
function by subtracting a polynomial function which mimics the weight function at those 
energies. In order to compensate for this subtraction, the same polynomial function is 
added again, but now its contribution is evaluated by using global duality on a circular 
contour in the complex plane, according to 

« 2 f S2 P / Ms) a 2 p (K(s) 

-L I U had (s)P n (s)ds U^ d (s)P n (s)d, 

ZTTZ J\a\=8i 2,1X1 J\s\=s 2 
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5 Application of the method to ajj ad 

The usual procedure to evaluate Eq. (J3|) is to use the experimental cross section up to 
some high momentum transfer and to calculate the remaining part from QCD using local 
duality. 

Indeed, in the energy region up to 1.2 GeV we use the e + e~ cross section data set from 
Refs. fS] and [TH]. The corresponding dispersion integral reads 



a 2 /■(1.2GeV)3 D BW^1#(*) 



p SlT z J Ami 1 J s 



Narrow resonances are added in explicit form. In Eq. (|13|) the missing resonances are 
parametrized by a Breit-Wigner fit i? BW (s) and added to the remaining data sets. 

In the subsequent interval from 1.2 to 3.0 GeV the relevant e + e~ data from Refs. |0J, 
combined with data sets from the BES Collaboration can be efficiently replaced by 
theoretical input from QCD. By subtracting a polynomial function from the weight function 
K(s)/s we reduce the experimental input in this energy interval. Since the polynomial 
function is an analytic function in the whole complex plane, the remaining contribution can 
be evaluated off the real axis on a circular contour in the complex plane using perturbative 
and nonperturbative QCD results, 



where rif = 3 is taken for this interval. 

At first glance the polynomial degree n seems to be arbitrary. However, if n is chosen 
too large, both the theoretical error from the strong coupling a s (Mz) and from unknown 
nonperturbative contributions arising from higher order condensates lead to larger uncer- 
tainties. 

Moreover, we have to ensure that the calculated value is stable with respect of the 
variations of the lower polynomial degrees n where the less known contributions from 
condensates are negligible. This point will help us find an optimal energy domain for our 
method. 

An optimal choice for the polynomial degree should guarantee a good approximation 
of the polynomial function to the weight function and thus effectively replace the experi- 
mental data by theoretical input while reducing the total uncertainty from this region to 
a minimum. 
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in units of 10 10 


n — 1 


n = 2 


n = 3 


n = 4 




78.52 


74.59 


71.02 


68.37 




-8.63 


-3.06 


-0.93 


-0.51 


K ad ) 2 


71.60 


71.49 


69.21 


67.08 


data contribution 


12.1% 


4.3% 


1.3% 


0.8% 



Table 1: Contributions to (a^ ad )2 for different fitting polynomial functions P n (s) with 
degree n. The purely experimental and theoretical contributions according to Eq. (|14|) 
are listed separately. Because the error estimate becomes worse for n = 3 (cf. Table |2J), 
polynomial degrees n = 1 and n = 2 are used to obtain a mean value and a methodical 
error estimate, (71.55 ± 0.13) x 10~ 10 . 

Above 3.0 GeV the e + e~ data in Ref. [011121 show low uncertainties and can be integrated 
directly, 

{aT)3= f 2 r^° [^>) +B - w] £M d , (15) 

P 371^ Js 2 =(3.0GeV) 2 L J S 

In the last step the high energy tail starting from 40 GeV can be calculated from theory 
using local duality. This is reasonable since there are no resonance contributions above 
40 GeV. In this region it is sufficient to use the lowest order contribution given by 

ta tf- w = » E ^f?! 1 ^)' (16) 

flavours f ±Z/I V 5 V b / 

6 Numerical results 

Using the recent e + e~ data from Refs. jH] and ^Hj we obtain for the dominant low energy 
component of a^ ad , 

(a**)! = (518.1 ± 4.8) x 10~ 10 . (17) 

We have added the statistical and point-to-point systematic errors in quadrature. In ad- 
dition to this experimental value we have to consider contributions from the u and 



8 



in units of 10 10 


n — 1 


n = 2 


n — 3 


n = 4 


uncert. due to a s {Mz) 


2.5 


5.1 


7.7 


10.5 


uncert. due to ((a s /7i)G 2 ) 


0.044 


0.196 


0.512 


1.043 


uncertainty of (a^ hc )2 


2.5 


5.1 


7.7 


10.6 



Table 2: Theoretical uncertainties of (cz^ ad )2 from QCD parameters and condensates. Be- 
cause the uncertainty is worse for higher polynomial degrees, the degrees n = 1 and n = 2 
are selected for the analysis. The error estimate is given by the average value 3.80 x 10 -10 . 
The uncertainty due to a s \(qq)\ 2 is less than 10~ 14 in all cases and is therefore omitted in 
the table. 

resonances for which we obtain the values (38.9 ± 1.4) x 10~ 10 and (40.4 ± 1.3) x 10 -10 , 
respectively, from an integration over the Breit-Wigner distribution function where the 
values for the parameters are taken from Ref. j20j . 

Applying our polynomial technique for the range between 1.2 and 3.0 GeV, we incor- 
porate QCD-expressions for the hadronic current-current correlator with corresponding 
uncertainties for the parameters occuring in these expressions. The errors from the masses 
of the light quarks u, d, and s can be neglected at energies above 1 GeV. Thus the theoret- 
ical uncertainty is dominated by the strong coupling a s (Mz) whose error originates from 
the uncertainty in the QCD scale = (380 ± 60) MeV. For the condensate terms we 
assign generous errors of 100%. Thus we have 

{—G 2 ) = (0.04 ± 0.04) GeV 4 , a s \(qq)\ 2 = (4 ± 4) x 10" 4 GeV 6 . (18) 

7T 

With increasing polynomial degrees the condensate error contributions grow fast. However, 
as can be clearly seen in Tab. ^ or in Fig. good approximations to the kernel of the dis- 
persion integral and consequently a very high reduction of the experimental data influence 
in the interval between 1.2 and 3.0 GeV can be achieved even by low degrees n = 1 and 
n = 2. This can be understood by the fact that one is far away from the singularity at 
s = and thus the weight function can be well approximated by low degree polynomials. 
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Figure 2: The QED kernel K(s)fs with fitting polynomial functions P n (s) of different 
degrees n. The polynomial functions were calculated by performing a least squares fit with 
no further constraints. 

We thus take the values (a^ ad ) 2 for n = 1 and n = 2 and calculate the algebraic mean. 
Therefore, we have to take into account a small methodical error due to the variation of 
the polynomial degree, 

(a™ eth ) 2 = ±0.06 x 10~ 10 . (19) 

Summing up the error contributions in quadrature and the algebraic mean uncertainty (for 
n — 1 and n = 2) from Table 121 one obtains 

( a J ad ) 2 = (71.55 ± 3.83) x 10~ 10 . (20) 

In the interval following 3.0 GeV up to the end of the e + e~ data set at 40 GeV as given 
in Ref. |Hj we obtain 

( a J ad ) 3 = (23.2 ± 0.2) x 10" 10 . (21) 

The contributions from narrow charmonium and bottonium resonances, (8.8 ±0.6) x 10~ 10 
and (0.11 ±0.01) x 10~ 10 , respectively, are taken by integrating over the narrow resonance 
distribution function where values for the parameters are again taken from Ref. [2~U] . 

10 
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Figure 3: Comparison of the l.h.s. and the r.h.s. of the sum rule Eq. ITU in the interval 
between 1.2 GeV and 3.0 GeV. Dotted horizontal line: value of integrating the l.h.s. using 
experimental data including error bars. The points give the values of the r.h.s. integration 
for various orders n of the polynomial approximation. Straight line interpolations between 
the points are for illustration only. The dashed lines indicate the error of our calculation. 

We finally obtain a tiny continuum contribution from the last interval, 

( a J ad ) 4 = — [°° T m n had (<0^U = 0.15 x 10~ 10 . (22) 



7T J (40 GeV) 2 S 

Summing up the contributions from the different energy intervals as shown in Table El 
we finally obtain 

a had = ( 701 3 ± 6 4 ) x 10 -io (23) 

for the hadronic contribution to a^. We have thereby assumed that errors with different 
origin are uncorrelated and have to be added quadratically in order to obtain the total 
error. 
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interval for ^/s 


contributions to a^ ad 


comments 


[0.28 GeV, 1.2 GeV] 


(518.1 ±4.8) x 10~ 10 


e + e~ cross section data 


u resonance 


(38.9 ± 1.4) x 10~ 10 


Breit-Wigner 


resonances 


(40.4 ± 1.3) x 10~ 10 


narrow resonances 


[1.2 GeV, 3.0 GeV] 


(71.6 ±3.8) x 10" 10 


polynomial method 


J/ip resonances 


(8.8 ±0.6) x 10~ 10 


narrow resonances 


[3.0 GeV, 40 GeV] 


(23.2 ±0.2) x 10~ 10 


e + e~ annihilation data 


T resonances 


(0.11 ±0.01) x 10- 10 


narrow resonances 


[40 GeV, oo] 


0.15 x 10- 10 


theory 


top quark contr. 


< 10~ 13 


theory 


hadronic contr. 


(701.3 ±6.4) x 10- 10 





Table 3: The different contributions to the hadronic part of the anomalous magnetic 
moment a^ ad of the muon. 

7 Higher order correction 

In the next-to-leading order (NLO) we have to consider three types of diagrams, those of 
type (2a) with an additional photon exchange, those of type (2b) with an electron loop 
inserted in one of the photon lines in Fig. ^ and finally the one of type (2c) with two 
correlator functions included in the photon line. For the former two the contribution reads 

or (NLO)4(^) 3 £^ W «M. (24) 

For numerical purposes it is convenient to represent the kernel functions K^ 2a \s) and 
(s) in terms of power series expansions in terms of m 2 / s [23] {m = rn^ = 105.6583568± 
0.0000052 MeV t 20j is the mass of the muon). One has 

K 1 s \ V 54 3 36 Km 2 ) J 

m 2 /8785 37tt 2 367 , / s \ 19 , 2 / s \\ 
+ v(ll^-^-216 ln UJ + l44 ln Km 2 )) (25) 
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m 4 / 13072841 883vr 2 10079 / s \ 141 
+ 7 2 " \ 432000 24CT ~ 3600 R W/ + ~80~ n ~ 



m 2 



*"w-4{(-n+i-fe)) 



7, 

m 2 / 55 7T 2 5, / s\ 5, /m 2 \ 1 2 /s\ l l2 /™ 2 \\ 
+ T(-48 + 18 + 9" 1 l^j + ^ ln M"6" 1 Ui + 6 ln UJJ <26) 

s 2 V 1800 3 3 \m) ) 10 \ m / / \ m f J \ m f J J J 

where for mj we used the electron mass, mj = 510.998902 ± 0.000021 keV [20j, since the 
contributions from the muon is already included in the contribution (2a). The contribution 
of the r is suppressed by m 2 /m 2 and numerically negligible. The kernels of the different 
orders in the electromagnetic coupling, K(s) on the one hand and K {2a \s) and K^ 2b \s) 
on the other hand, are compared in Fig. 0] For illustration, we like to sum the results from 
the evaluation of Eq. (|2*4*j) which was performed in a completely analogous way as in the 
previous chapter. The final result (-21.07 ± 0.21) x 10~ 10 in case of K^ a \s) (cf. Table ij) 
and (10.78 ± 0.08) x 10 -10 in case of K^ 2h \s) (Table 0) are in good agreement with the 
literature (221- Finally, for the small contribution from diagrams of type (2c) with two 
correlator functions included into the photon line we take the value (0.27 ± 0.01) x 10~ 10 
from Ref. [21]. We have checked that to the required accuracy a recalculation using the 
method presented in this paper is not necessary. The total hadronic NLO contribution then 
sums up to a value of (—10.3 ±0.2) x 10~ 10 , as compared to the value (—10.1 ±0.6) x 10 -10 
given in Ref. [23J. Our central value is 3% smaller than that given in Ref. [22] and the error 
is reduced from 0.6 x 10~ 10 to 0.2 x 10~ 10 . Finally, summing up leading and next-to-leading 
order contributions together with the afore-mentioned QED and weak contributions and 
the so-called light-by-light contribution a]} ad (LBL) = (8.6 ± 3.5) x 10~ 10 
obtain a M = (11 659 185.6 ± 6.4 had ± 3.5 LBL ± 0.4 QED+EW ) x lO" 10 . 
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interval for yHs 


contributions to a^ d 


comments 


[0.28 GeV, 1.2 GeV] 


(-14.04 ±0.12) x 10- 10 


e + e~ cross section data 


uj resonance 


(-1.11 ±0.04) x 10- 10 


Breit-Wigner 


resonances 


(-1.29 ±0.04) x 10- 10 


narrow resonances 


[1.2 GeV, 3.0 GeV] 


(-2.84 ±0.14) x 10~ 10 


polynomial method 


J/ijj resonances 


(-0.44 ±0.03) x 10~ 10 


narrow resonances 


[3.0 GeV, 40 GeV] 


(-1.35 ±0.09) x 10~ 10 


e + e~ annihilation data 


T resonances 


< io~ 13 


narrow resonances 


[40 GeV, oo] 


< io- 12 


theory 


top quark contr. 


< 10~ 14 


theory 


hadronic contr. 


(-21.07 ±0.21) x 10~ 10 





Table 4: The different contributions to the hadronic part of the anomalous magnetic 
moment a^ ad of the muon for the kernel K^ 2a \ 



interval for y/s 


contributions to a*f d 


comments 


[0.28 GeV, 1.2 GeV] 


(7.98 ±0.07) x IO" 10 


e + e~ cross section data 


u resonance 


(0.59 ±0.02) x IO" 10 


Breit-Wigner 


resonances 


(0.62 ±0.02) x 10~ 10 


narrow resonances 


[1.2 GeV, 3.0 GeV] 


(1.09 ±0.01) x 10~ 10 


polynomial method 


J/ip resonances 


(0.14 ±0.01) x 10~ 10 


narrow resonances 


[3.0 GeV, 40 GeV] 


(0.364 ±0.003) x IO" 10 


e + e~ annihilation data 


T resonances 


< IO" 13 


narrow resonances 


[40 GeV, oo] 


< IO" 12 


theory 


top quark contr. 


< io- 14 


theory 


hadronic contr. 


(10.78 ±0.08) x 10~ 10 





Table 5: The different contributions to the hadronic part of the anomalous magnetic 
moment a^ ad of the muon for the kernel K^ 2b \ 
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0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

s [GeV 2 ] 

Figure 4: Comparison of the LO and NLO kernels K(s) and K^ 2a ' 2b \s) normalized at 
s = 4m 2 . The functional behaviour of the three weight functions is quite similar. 

8 Conclusion 

We have presented an alternative determination of the hadronic contribution to the anoma- 
lous magnetic moment of the muon where we have made use of theoretical QCD results to 
reduce the influence of the poor experimental data in the range between 1.2 and 3.0 GeV. 
Our analysis includes the leading and next-to-leading contribution. 

Using the polynomial method in the range between 1.2 and 3.0 GeV, we could suppress 
the influence of experimental data effectively and thereby reduce the error on the determi- 
nation of a^ ad . Note that using this method we were able to use a value for the lower limit 
of this range lower than usually used for pure QCD methods. 

Using only e + e~ data and QCD as input, we obtain a result which is 1.6a away from 
the measured world average, stating that the deviation between theory and experiment 
may be smaller than commonly suggested. Other recent results to compare with are 
(11659180.9±7.2 had ±3.5 L BL±0.4 Q ED+Ew)xl0~ 10 |2ZJ and (11 659 166.9±7.4) x 10~ 10 0. 
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